x = (20:2000:125010)';
y = lognpdf(x,log(20000),1.0);
plot(x,y)
